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Abstract 

Contact  theorems  for  rough  surfaces  are  discussed.  Simple  relations 
for  the  average  contact  density  are  obtained  for  neutral  and  charged  walls. 
When  the  walls  are  not  planar  there  are  new  contributions  proportional  to 
the  held  gradients  near  the  charged  wall. 


1  Introduction 


Contact  theorems  are  very  useful  to  estimate  the  accuracy  of  integral 
equations  and  computer  simulations  for  inhomogeneous  systems  near  neutral 
or  charged  walls. 

We  have  derived  in  the  past  sum  rules  for  liquids  near  charged  plemar 
surfaces  [1,  2]  in  which  the  total  contact  density  was  given  as  a  function  of 
the  charge  on  the  wall  and  the  bulk  pressure  of  the  fluid. 

In  recent  times,  and  due  to  the  rapid  increase  in  computer  performance, 
simulations  of  more  realistic  systems  are  possible.  Non  spherical  solute 
molecules,  and  metal  surfaces  with  crystal  structures  are  accessible  to  sim¬ 
ulations  of  interfaces  with  discrete  (and  also  non  spherical)  solvent  (water) 
molecules.  In  these  simulations  it  is  very  important  to  use  a  potential  with 
a  correct  equation  of  state  (pressure).  In  general,  in  the  modeb  of  water  the 
interaction  is  either  adjusted  to  the  experimental  value  of  the  density,  and 
there  is  an  error  in  the  pressure,  or  (in  a  few  cases)  the  interaction  is  ad¬ 
justed  to  the  pressure,  and  there  is  a  generally  slight  error  in  the  density.  For 
the  simulation  of  interfaces,  we  think  that  the  latter  situation  is  preferable. 
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The  contact  sum  rules  are  based  on  force  balance  considerations,  between 
the  bulk  pressure  and  the  forces  exerted  by  the  molecules  on  the  container 
walls.  Although  a  formal  derivation  could  be  obtained  by  integration  jf  the 
Born-Green-Yvon  equation,  it  can  be  shown  that  all  pair  (or  higher  multi- 
plet)  interactions  of  molecules  in  the  fluid  phase  cancel  and  therefore  only 
singlet  densities  count.  Our  results  are  extensions  of  previous  work  [3,  4,  5]. 

There  are  two  kinds  of  dynamic  relations  that  are  interesting  in  electro¬ 
chemistry:  those  for  the  normal  pressure  and  those  for  the  stress  and  strain 
along  the  interface,  that  have  been  measured  experimentally  [6].  We  discuss 
in  this  paper  the  first  ones  for  surfaces  that  are  planar  on  the  average,  but 
not  necessarily  locally.  This  includes  rough  surfaces,  but  more  interestingly 
interfaces  between  a  single  crystals  and  solutions. 

The  most  general  definition  of  a  planar  interface  is  obtained  by  requir¬ 
ing  that  the  average  force  acting  on  an  individual  molecule  of  the  fluid,  of 
species  i  =  1 . . .  m  is  integrated  to  a  constant  for  arbitrarily  large  surfaces. 
The  surface  defined  by  zi{x,y]A)  represents  the  contact  surface,  and  is  the 
solution  of  the  equation 

Wi{x,y,z)^A  (1) 

where  toj(r)  =  Wi(x,y,z)  is  the  interaction  of  the  wall  with  molecule  i  for 
given  x,y,  and  finite  A.  Then 

^  j  dxdyzi{x,y]A)  =  Bi  (2) 

and  our  condition  simply  is  that  eq.(2)  is  satisfied  for  every  component  of 
the  mixture.  This  is  a  planar  surface  on  the  average. 

For  hard  core  potentials  we  clearly  choose  zi{x,z)  to  be  the  hard  core 
surface,  and  the  requirement  is  that  the  average  z\  over  the  entire,  infinite 
surface  is  a  constant. 


2  Neutral  systems 


Consider  the  total  energy  of  the  interface 

i 


(3) 
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where  lOj(l)  =  Wi{fi-,  {6},  Jli)  is  the  total  interaction  potential  of  the  wall  and 
a  particle  of  species  i  situated  at  position  fi  with  orientation  flj  .  (fii  are  the 
Euler  angles  that  define  the  molecular  orientation).  {6}  is  a  set  of  parameters 
which  define  the  position  of  the  surface  :  if  we  produce  an  arbitr<iry  but  small 
deformation  of  the  surface,  we  get  from  eq.(3)  the  Taylor  expansion 


SW  =  J 


We  are  interested  in  a  particular  deformation  of  the  interface:  a  uniform 
displacement  in  the  z  direction:  Then 


P  = 


Po  =  - 


1  dW 
S  dz 


^  J^dz,dy„m,  Jd^,'£Ml)^ 


(5) 


which  expresses  the  force  balance  equation  in  the  z-direction. 


Let  us  consider  first  the  case  of  hard  walls  near  a  fluid  with  hard  con¬ 
vex  molecules.  The  simplest  case  is  that  of  a  planar  wall.  The  hard  core 
interaction  is  best  expressed  by  the  force  relation 

^^  =  -»:sr<[Fi(f,,f!,)-l)  (6) 

where  the  solution  of  F«(ri,ni)  =  1  yields  the  position  of  contact  of  the 
molecule  at  the  wall  [8],  when  the  orientation  is  given  by  the  Euler  angles 
<^i)-  Then  we  obtain  Zi(f2i)  as  the  solution  of 

Fi(zi,Sli)  =  l  (7) 

For  a  hard  ellipsoid,  of  main  axes  A,  B,  C  for  instance  we  get  [8] 

=  ^l/{cos®  ^,/C*  +sm^5i[sin*  dilA?  +  cos*^/B*]}  (8) 

where  Rij  are  the  matrix  elements  of  Euler’s  rotation  matrix. 
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Replacing  into  eq.(5),  we  get 

P  =  p^  =  kBT^j  (9) 

where  21(^1)  is  given  by  eq.(8)  When  A  =  B  =  C  we  get  from  eq.(8) 

zi  =  A  =  ct/2  (10) 

and  we  recover  the  well  known  relation 

P=P„  =  kTY,M2)  (11) 

i 

When  the  walls  are  not  flat,  then  we  have  to  discuss  the  normal  com¬ 
ponents  of  the  pressure  tensor,  and  the  tangential  components  of  the  stress 
tensor.  Consider  only  crystal  surfaces,  which  aire  more  interesting  experimen¬ 
tally.  The  equation  for  the  contact  distance  is  more  compHcated,  since  we 
have  to  solve  a  rather  more  complex  equation  involving  convex  and  concave 
regions 

Fiixiziyi;ni) -1  =  0  (12) 

which  is  solved  to  yield 

zi(l)  =  z(ii,yi,ni)  (13) 

To  compute  the  force  we  need  to  know  the  equation  for  the  normal  to  the 
tangent  plane  more  specifically 

cos/3,(xiyi)  =  n,.e,  (14) 

where  Uj  is  the  direction  of  the  tangent  plane  of  the  molecide.  As  an  example 
consider  the  triangular  lattice  of  lattice  constant  d  and  height  h:  The  simplest 
form  of  the  surface  is  [7] 

z{x,y)  =  {h/ 2)  [cos  a  +  cos  b]  (15) 


where 


a  =  —{x  -yiVS)  and  6=— (2y/\/3) 
a  a 


4 


4 


The  components  of  the  noimal  vector  n,  =  n,{x,  y)  are  easily  computed  from 
dzjdx  and  dzidy: 


cos^,  = 


1 


1  +  {7r/i/d)^[sin^  a  +  (l/3){sina  —  2  sin  6)^] 
The  normal  pressure  sum  rule  is  then 


Po  =  fcflT  ^2  ^<(con) 


(16) 


(17) 


where 


Pi{con)  =  ^  J  dxidyi  j  -^pi[zt{l)]cos /S,{xi,yi)  cose„  (18) 

where  now  $„  is  the  angle  between  the  tangent  plane  at  contact  and  the 
direction  of  the  line  joining  the  contact  point  and  the  center  (also  center  of 
mass)  of  molecule  i. 


3  Charged  Systems 

The  interaction  uj,(r)  is  now  the  sum  of  two  contributions:  a  Coulomb 
or  electrostatic  term  plus  a  non  electrostatic,  covalent  term  such  as  discussed 
above. 

u)i(r)  =  u;-®*'(r)  +  t»'*(r)  (19) 

Consider  different  situations: 

1.  The  uniformly  charged  electrode  facing  a  primitive  model  (continuum 
dielectric)  electrolyte.  In  this  case 

Wi{r)  =  w^{r)  +  Wi’ir)  (20) 

The  electrostatic  contribution  is  in  this  cjise 

u;'*(z)  =  -eiEozjl  (21) 


5 


where  is  the  chaxge  of  i,  £^o  is  the  bare  electric  field.  Using  the 
electroneutrality  relation 

where  qs  is  the  surface  charge  density  of  the  elctrode.  We  immediately 
get 


P  =  Po  -  ~[E,Y  (23) 

where  Po  is  given  by  Eq.(17). 

It  expresses  the  fact  that  the  total  pressure  must  be  equal  to  the  ki¬ 
netic  term  due  to  the  collisions  of  the  molecules  at  the  wall  minus  the 
attractive  electrostatic  contribution  of  a  planar  capacitor  with  charge 
density  qs  =  ^. 


2.  The  flat  electrode  facing  a  non  primitive  model  electrolyte.  If  the  sol¬ 
vent  consists  of  neutral  molecules  with  a  dipole  (  or  higher  multipoles 
),  there  will  be  no  net  force  since  the  dipoles  interact  with  the  gradient 
of  the  applied  field  VEo,  which  in  this  case  is  zero.  Therefore,  only  the 
hard  repulsive  interactions  count,  and  we  get  now 

P  =  Po{solute)  -f  pQ{solvent)  —  (24) 

OTT 

where  Po(so/ute),  Po(so/uent)  are  the  bulk  pressure  terms  given  by 
Eq.(17).  Notice  that  now  the  dielectric  constant  has  disappeared  from 
this  relation.  This  means  that  the  electrostatic  contribution  in  a  sol¬ 
vent  of  high  dielectric  constant  like  water,  is  now  much  smaller  than 
in  the  primitive  model,  and  that  the  covalent  (heird  core)  solvent  term 
Po  plays  a  much  larger  role  in  the  makeup  of  the  concentration  profile 
near  the  electrode. 

3.  The  rough  or  structured  electrode  near  a  primitive  electrolyte.  The 
situation  is  now  more  complicated  since  the  charge  distribution  at  the 
electrode  surface  will  not  be  uniform,  amd  therefore,  both  contact  den¬ 
sities  (for  hard  surfaces  as  well  as  for  soft  surfaces)  will  be  functions 
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of  z  as  well  as  the  position  on  the  surface  x,y.  There  will  be  a  simple 
relation  only  for  average  quantities  such  as  the  average  contact  density 
near  a  hard  plane,  Eq(18) 

It  is  clear  that  the  electrostatic  forces  along  the  surface  are  of  vanishing 
magnitude  for  a  random  rough  surface,  or  zero  for  a  periodic  crystal 
surface.  We  assume  that  the  surface  has  a  fixed  charge  distribution, 
and  since  we  are  neglecting  images,  it  is  not  metallic. 


Asume  now  a  capacitor  in  vacuum,  in  which  one  plate  is  our  surface  S, 
with  its  charge  distribution,  and  the  second  plate  is  very  far,  and  has  a 
uniform  charge  density  —qs,  such  that  there  is  no  field  sufficiently  far 
behind  the  first  plate  of  our  capacitor.  From  Gauss’s  theorem,  and  since 
the  transversal  components  of  the  field  E2(l),  E°(l)  are  as3Tnptotically 
zero  at  the  sides  of  our  capacitor 

<  E"  >s=  (1/5)  J  dxtdy^E°{l)  =  {iw/e)qs  (25) 

Since  the  far  away  electrode  creates  a  field  —(1/2)  <  E°  >s,  which  does 
not  contribute  to  the  pressure  on  the  first  electrode,  the  electrostatic 
contribution  to  the  pressure  P*  is  given  by 

P.  =  l/S  fdlf;  [£?(!)  -  (1/2)  <  £?  >sj  = 


/ni  M  m 

dlY,“PMW2)  <  >s+l/S  /  <il5^ei«(l)lf:.”(l)-<£j>s] 
.=1  ''  1=1 

(26) 

Using  charge  electroneutrality,  we  get 


n  =  <  £?  >1  +1/5  jdl'^ .,ft.(l)l£j(l)-  <  E°  >s]  (27) 

where  the  first  term  is  the  Maxwell  stress  tensor  contribution  present 
in  the  flat  electrode  case,  while  the  second  term  is  a  new  one  arising 
from  fluctuations  in  the  bare  field.  The  integration  should  cover  only 
a  narrow  region  near  the  electrode  surface. 
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The  other  single  particle  term  contjiinig  the  short  range  interactions 
between  the  molecules  and  ions  and  the  wall,  yields 


<  >s=  US  J^dx^dy,pi(l)s7,  u^r(l)  (28) 


where 


<  PiiU — ^ - ->s=l/S  I  dxidyi  I  dzipi 


Putting  it  all  together  )delds  the  general  contact  theorem  for  a  planar 
on  the  average,  but  not  necessarily  smooth,  surface 


P  =  Pc 


>S 


+1/5  /  Jlf^e,ft(l)l£j(l)-  <  Ef  >s] 
•'  .=  1 


-  E  <  '■‘(i) 

i=i 


dzi 


>s 


(30) 


This  theorem  is  a  generalization  of  the  previously  derived  contact  theo¬ 
rems  to  the  realistic  case  of  non  smooth  electrode  surfaces.  It  contains 
the  previous  results  as  partictilar  cases. 


For  a  surface  vdth  an  array  of  sticky  adsorption  sites,  such  as  in  the 
case  of  the  sticky  site  model,  (SSM)  model  discussed  elsewhere  [9]  ,  the 
adsorption  potential  has  the  form 


with 


-/3«i«(rT  ^  ^  X„{R)6{z) 

(31) 

M 

=  2^  KS{R  —  niSi  —  0202) 

(32) 

T»i  ,nj 


Here  R  =  x,y  is  the  position  at  the  electrode  surface,  and  z  the  distance 
to  the  contact  plane,  which  is  at  a  distance  (t/2  irom  the  electrode.  In 
32,  n\ ,  n2  are  integers,  there  are  M  sites  on  the  electrode  of  area  5  ,  and 
Si ,  02  aie  the  lattice  vectors  of  the  adsorption  sites  on  the  surface.  The 
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parameter  Ag  represents  the  fugadty  of  an  adsorbed  atom  of  species  a. 
Define  now  the  regular  part  of  the  density  function 


ViW  =  (33) 

Replacing  into  the  general  contact  theorem  30  gives  [5] 

P  =  IgT  ft(<T,/2)  -  6/8»(£;(0)P  +  ^  ^  >  P‘  (34) 

1=1  1=1  ' 

This  theorem  has  been  verified  recently  by  [4],  for  the  exactly  solved 
model  of  a  one  component  plasma  in  two  dimensions. 

4.  The  rough  electrode  near  a  non  primitive  {for  example  with  a  solvent 
of  dipolar  hard  spheres)  electrolyte.  Now  we  have  to  include  the  effect 
of  electric  field  gradients,  which  are  not  zero  near  the  electrode.  The 
total  electrostatic  force  is  [10] 

-  =  eJo  +  MVEo)  +  (1/6)?,  :  (VV^o)  (35) 


where  is  the  dipole  moment  of  a,  qa  its  quadrupole  moment,  and 
so  on.  We  remark  that  now  the  solvent  singlet  peirticle  density  Pa{l)  is 
not  only  a  function  of  fi,  but  also  of  the  orientation  of  the  molecules 
with  respect  to  the  electrode,  which  in  the  case  of  the  lineeir  dipoles  is 
given  by  <^i.  Therefore  we  expand 

/>.(!)  =  E =  P°.o  +  +  -  (36) 

l,m  m 


The  dipole  contribution  to  the  pressure  is 


4ir^ 


i  jf”  <(z.(l/5)  J  (37) 


which  after  a  short  calculation  leads  to 
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P  =  pQ{solute)  +  PQ{solvent)  —  ^  <  ^*(0) 

oTT 


47r/i.  .  .1  9E\^{fx) 


+1/S  J  <  £»(o)  >i-i^  <  '£p'.. 

«=l  m 


dzi 


>. 


where  we  have  used 


(38) 


=  (39) 

m 

where  are  the  polar  components  of  the  unit  vector. 

The  last  term  of  Eq.(38)  corresponds  to  an  electrostriction  effect  which 
vanishes  for  uniform  external  field  Eq. 

5.  The  rough  electrode  near  a  non  primitive  electrol3d;e:  Most  of  the  mod¬ 
els  of  water  employed  in  the  computer  simulations  consist  of  neutral 
molecules  with  embedded  point  charges. 

The  sum  over  the  charges  q,,  in  each  molecule  is  indicated  by  the  index 
I/,  and  is  zero  for  each  molecule.  Each  of  these  charges  is  located  at 
the  position  relative  to  a  molecular  reference  frame.  From  Eq.(30) 
we  get 


P  =  Po{solute)  +  Po{solvent)  —  ~  <  -£*(0)  >* 

oTT 

+1/S  I  <ilf;e,A(l)[f;(l)-<B«(0)>]+ <;,,(!)  ,„£5(f,+6„)  >s 

1  =  1  v 

(40) 

where  p*(l)  is  now  the  orientation  dependent  density  function  of  the 
solvent  with  embedded  charges  at  contact. 

These  relation  point  to  the  importance  of  using  a  model  potential  for  bquid 
water  that  has  the  correct  equation  of  state  (pressure)  rather  than  the  correct 
bulk  density  (  which  may  be  off  may  be  by  a  few  percent)  when  computing 
density  profiles  near  planar  or  rough  electrodes. 
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